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disappointing.  However,  when  the  elc"  ated  duct  is  allowed  to  rise  slowly  with  increasing  range  the 
RK)  results  closely  match  the  observed  signal  levels,  indicating  (hat  weakly  alicnuaicd  l.atth 
detached,  or  whispenng  gallery,  modes  were  the  dominant  propagation  mechanism 


INTRODUCTION 

Papperi  and  Goodhart  (1977]  presented  a  case 
study  of  propagation  at  3087.7  MHz  in  the  presence 
of  an  elevated  duct.  In  this  case  study  the  transmit¬ 
ter  and  receiver  heights  were  21  and  914  m  above 
sea  level,  and  measurements  were  made  to  about 
700  km.  The  transmitter  utilized  a  horizontally 
polarized  standard  gain  horn  antenna  with  a  beam 
width  of  32°  on  the  coast  of  Point  Loma  in  the  city 
of  San  Diego.  The  receiver  utilized  an  open-ended 
waveguide  feed  antenna  with  a  beam  width  of  85°. 
The  receiver  and  data  collection  systems  were 
flown  in  an  aircraft  on  an  entirely  overwater  path 
westward  from  the  transmitter  site.  The  case  se¬ 
lected  was  May  28,  1974,  and  the  measurements 
were  made  on  an  outbound  flight  from  1853  to  2140 
UT  (1153-1440  PDT).  Figure  1  shows  modified 
refractivity  profiles  measured  by  radiosonde  at  the 
Point  Loma  transmitter  site  at  1909  and  2247  UT. 
Modified  refractivity  in  M  units  is  given  by 

M  =  10*(/i  -\+zJa) 

where  n  is  the  refractive  index,  z  is  height  above  sea 
level,  and  a  is  the  Earth's  radius.  Both  profiles 
indicate  that  an  elevated  duct  existed  from  approx¬ 
imately  350  to  700  m  above  sea  level.  Figure  I  also 
shows  the  transmitter  and  receiver  heights,  which 
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vycre  substantially  belotv  and  above  the  elevaied 
duct. 

f  igure  2.  reproduced  from  Figure  1 7  of  Papperi 
and  Goodhari  11977]  (with  a  legend  conection). 
compares  path  loss  in  decibels  versus  range  for  ihe 
measured  data  with  results  calculated  using  a 
waveguide  model  and  a  trilinear  refractivity  profile 
based  on  the  measured  profiles  of  Figure  1.  Refer¬ 
ence  curves  for  free-spacc  and  standard  atmosphere 
losses  are  also  included  in  Figure  2.  The  waveguide 
results  match  the  mea.sured  data  quilc  we))  at  near- 
horizon  ranges  but  ovcrprcdict  the  loss  by  more 
than  60  dB  at  the  greater  range,s.  Pappen  and 
Goodhart  observed  that  the  measured  data  arc 
indicative  of  energy  being  transported  by  weakly 
attenuated  modes  and  speculated  that  mode  conver¬ 
sion  proces.scs  from  a  range-dependent  environ 
ment  might  be  important  in  this  case  study.  They 
also  suggested  that  a  strong  low-lcvcl  evaporation 
duct  might  account  for  the  discrepancies  between 
the  modeled  and  measured  loss  values. 

Hitney  et  al.  [1978]  reexamined  this  case  study  in 
terms  of  evaporation  duct  propagation.  Since  the 
evaporation  duct  was  not  expected  to  affect  the 
relatively  low  microwave  frequency  of  the  experi¬ 
ment,  accurate  meteorological  measurements  near 
the  sea  surface  were  not  made  in  the  1974  experi¬ 
ment.  Therefore  Hitney  et  al.  used  hypothetical 
multilevel  refractivity  profiles  reprc.senting  moder¬ 
ate,  strong,  and  very  strong  evaporation  ducts  with 
duct  heights  of  12,  23,  and  35  m.  respectively,  in 
their  study.  These  profiles  were  used  as  input  to  a 
waveguide  model  similar  to  the  one  described  by 
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Fig,  1 .  Modified  refractivny  profiles  measured  h\  radiosonde 
on  May  28,  1974,  showing  [he  eievaled  duel  and  the  iransmiUcr 
and  receiver  heights 


sirriahh  fjmii!  diouh  (ii  !iu  rm  ..‘-.iin  jiiiTi! •,  .iruJ  ihis 
Jlisl  s*li.’lii,!!ll  is  sjUlic  I. Ill  Hi  llu:  S,ili  .tic.i 

I'oi  csaiiipic,  lusiiHis.il  li.ti.i  liuin  Use  s-iipmcci  s 
(t'lcaCU's'  cflccls  pu-iiis  luni  sSsU-rn  j /’a.’ii  r  s.o;  t'l 
al  ,  !9‘X)1  shoss  a  '‘'  Hi  cvapoi.iiuiii  slus!  is  c\ps*stcsi 
tar  less  than  T  id  ihc  iiiiic  un  .ui  annua!  hasp,  m  itic 
San  Diego  aica  l  lius  il  iiuist  be  sonsluded  ihai  she 
evaporation  duct  is  a  possible  hui  highb.  unlikeK  . 
explanation  for  the  higfici -than  expes  led  signal  lev 
els  measured 

RADIO  HMYSK  Al  OPI  K  S  MODI  I 

I'he  radio  physical  optics  (RPOi  propagation 
mrxicl  is  a  range-dependent  hvhrid  model  (hat  com 
bines  ray  optics  (RO)  and  parabolic  equation 
methods  to  produce  a  relatively  fasi-runmng  com 
puter  program.  The  !^T  mode!  is  based  ors  rnethtHJs 


Pappert  and  Goodharl  (19771,  and  calculated  re¬ 
sults  v'ere  compared  to  the  experimental  measure¬ 
ments.  Figure  3,  reproduced  from  Figure  3  of  //ii- 
ney  et  al.  [1978],  shows  the  results  of  these 
comparisons.  Figure  3  shows  that  evaporation 
ducts  can  affect  signal  levels  far  above  the  surface. 
However,  only  the  very  strong  duct  shows  a  rea- 
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Fig.  2.  Path  loss  versus  range  for  May  28.  1974,  saowing 
experimental  measurements  (solid  line)  and  calculations  from 
waveguide  methods  (dashed  line),  repriKliiced  from  Pappert  and 
Goodhan  11977). 


Fig.  Path  loss  versus  range  for  Mav  28.  1974,  showing 
experimental  measurements  (solid  lines)  and  moderate,  strong, 
and  ver>'  strong  evaporation  duct  models  (dashed  Imcs)  from 
waveguide  methods,  reproduced  from  Hirnev  n  al  [1978], 
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Fig  4  The  four  radio  physical  optics  (RK))  model  regions 


described  by  Tapper!  [1977]  and  Dockery  [1988]. 
One  disadvantage  of  PE  models  applied  to  micro- 
wave  problems  is  that  inclusion  of  elevation  angles 
above  a  fraction  of  a  degree  can  result  in  very  long 
computer  execution  times  and  large  memory  re¬ 
quirements.  RPO  overcomes  this  disadvantage  by 
using  a  variety  of  RO  methods  above  a  small 
limiting  elevation  angle  to  compute  propagation 
loss,  thus  keeping  the  larger  angles  out  of  the  PE 
solution.  In  the  RO  methods,  full  account  is  given  to 
focusing  or  defocusing  along  both  direct  and  re¬ 
flected  ray  paths  and  to  the  integrated  opti  'al  path 
length  difference  between  the  two  ray  paths  to  give 
precise  phase  difference  and  hence  accurate  coher¬ 
ent  sums  of  the  RO  components.  To  further  en¬ 
hance  the  speed  of  computation,  the  fast  Fourier 
transform  size  of  the  PE  model  is  variable  but 
limited  to  a  maximum  of  1024  points.  This  implies  a 
maximum  receiver  height  for  the  PE  model  that 
depends  on  frequency,  transmitter  height,  and  re- 
fractivity  profiles.  For  typical  applications  at  3  GHz 
the  maximum  PE  receiver  height  is  about  2500  m. 
For  points  above  this  height,  but  beyond  the  RO 
region,  an  extended  optics  model  has  been  devel¬ 
oped  that  uses  ray  optics  methods  that  are  initial¬ 
ized  by  the  PE  solution.  Finally,  for  elevation 
angles  above  5°  or  ranges  less  than  2.5  km,  a 
flat-Earth  model  is  used  that  ignores  all  refractive  or 
curved-Earth  effects.  Figure  4  summarizes  the  rel¬ 
ative  location  of  the  four  RPO  regions. 


There  arc  three  primary  advantages  of  the  hybnd 
RPO  model,  as  discussed  by  Hiiney  [  1992].  First,  it 
has  proven  to  be  from  25  to  100  times  faster  than 
pure  split-step  PE  methods  lor  many  practical 
cases,  with  no  loss  in  accuracy.  .Second,  RPO  can 
be  applied  to  very  wide  elevation  angles,  whereas 
pure  PE  methods  arc  generally  restricted  to  small  or 
moderate  angles  due  to  practical  limitations  on 
memory  size.  The  third  advantage  is  that  the  total 
memory  requirement  of  RPO  is  quite  modest,  need¬ 
ing  only  160K  on  IBM/PC -compatible  computers. 

COMPARISONS  OF  RK)  TO  MEASUREMENTS 

The  RPO  model  was  applied  to  the  case  study  of 
May  28,  1974,  using  a  combination  of  the  measured 
profiles  shown  in  Figure  1  and  assuming  horizontal 
homogeneity.  The  exact  refractivity  profile  used  by 
RPO  is  given  in  Table  1  in  terms  of  M  units  versus 
height.  The  results  are  presented  in  Figure  5  as 


TABLE  I.  Modified  Refraciivity  Vcr.sus  Height  Profile 
Based  on  a  Combination  of  the  Two  Measured  Profiles  of 
Figure  1 
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homogeneous  elevated  duct  environment  of  May  28.  1974  (solid 
line),  compared  to  experimental  measurements  (dotted  line) 
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propagation  loss  versus  range  at  914  m  compared  to 
the  measured  radio  data.  In  this  paper,  path  loss 
and  propagation  loss  are  equivalent,  both  being  the 
ratio,  in  decibels,  of  the  transmitted  to  received 
power  assuming  loss-free  isotropic  antennas.  The 
RPO  results  were  calculated  assuming  an  omnidi¬ 
rectional  antenna  for  consistency  with  the 
waveguide  results,  but  there  would  be  virtually  no 
difference  in  loss  at  the  greater  ranges  if  the  actual 
transmitter  antenna  pattern  had  been  used.  Refer¬ 
ence  lines  are  also  included  in  Figure  S  for  frcc- 
space  and  receiver-threshold  loss  values.  These 
RPO  results  show  a  strong  similarity  to  the  results 
of  Pappert  and  Goodhart  presented  in  Figure  2, 
which  indicates  that  the  RPO  and  waveguide  mod¬ 
els  are  equivalent  for  the  homogeneous  case.  Un¬ 
fortunately,  the  exact  profile  used  by  Pappert  and 
Goodhart  could  not  be  determined.  Hitney  [1992] 
has  shown  that  RPO  and  waveguide  methods  give 
identical  results  when  applied  to  identical  profiles. 

Since  no  refractivity  measurements  were  made 
away  from  the  transmitter  site  in  the  1974  experi¬ 
ment,  there  is  no  way  to  ascertain  how  the  refrac¬ 
tivity  profile  may  have  changed  along  the  propaga¬ 
tion  path.  However,  Figure  19  of  Hitney  et  al. 
[1985]  shows  a  duct  similar  to  those  of  Figure  1  that 
rose  with  a  slope  of  1:800  (that  is,  1  m  in  height  for 
every  800  m  in  range)  on  an  overwater  radial  away 
from  Point  Loma.  Also,  Figure  4.24  of  Ken  [1951] 
shows  a  layer  with  a  slope  of  about  1:1000,  and 
Figure  3  of  Guinard  et  al.  [1964]  shows  a  layer  with 
a  slope  of  about  1:1700,  both  of  which  were  ob¬ 


served  in  the  San  Diego  area.  A  discussion  of  the 
meteorological  conditions  that  result  in  these  rising 
ducts  is  given  by  Guinard  et  al.  A  few  cases  of  range 
dependency  were  therefore  modeled  by  RI  O  with 
constant  slopes  up  to  1:800  using  the  same  refrac¬ 
tivity  profile  at  the  transmitter  as  was  used  for 
Figure  5.  The  best  fit  of  these  range-dependent 
cases  was  one  with  a  slope  of  1.1500,  the  results 
from  which  are  shown  in  Figure  6.  This  case  used 
the  profile  from  Tabic  I  at  range  zero  and  the  profile 
from  Table  2  at  750  km.  Note  the  modeled  loss 
levels  at  the  far  ranges  are  in  good  agreement  with 
the  measured  data,  even  though  at  the  near  over- 
the-horizon  ranges  they  are  not  much  better  than 
the  homogeneous  case.  In  this  case,  allowing  the 
duct  to  rise  very  slightly  has  increased  the  signal 
levels  by  about  50  dB  compared  to  the  homoge¬ 
neous  case. 

The  reason  why  the  range-dependent  modeled 
signals  are  much  higher  than  the  homogeneous  case 
is  apparent  in  the  coverage  diagram  of  Figure  7, 
which  represents  the  same  case  as  Figure  6.  Figure 


TABLE  2.  Modified  Refractivity  Versus  Height  Profile 
Assumed  al  750  km  for  the  Rising  Duct  Case  With  .Slope  of 
1  1500 
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Fig.  7.  Coverage  diagram  from  RPO  for  the  elevated  duct  environment  of  May  78,  1974.  vmh  duel  rising  at  a 

slope  of  l  l-iOO. 


7  shows  propagation  loss  in  lO-dB  increments  on  a 
gray  scale  on  a  height-versus-range  display.  It  is 
apparent  that  an  aircraft  flying  level  at  914  m  would 
intersect  the  duct  at  ranges  greater  than  400  km, 
where  the  received  signals  would  be  dominated  by 
the  Earth-detached  modes,  or  whispering  gallery 
effects,  associated  with  the  elevated  duct.  The 
Earth-detached  modes  in  the  elevated  duct  are 
clearly  excited  by  the  source  that  is  well  below  the 
duct  and  propagate  with  very  low  attenuation  rates. 

A  few  other  range-dependent  refractivity  envi¬ 
ronments  were  modeled  by  RPO  to  see  if  the 
measured  signal  levels  between  about  150  and  300 
km  could  be  better  matched.  The  best  match  was 
found  using  a  duct  with  an  upward  slope  of  1:750 
between  0  and  250  km  and  no  slope  beyond  250  km. 
Again,  the  refractivity  profile  at  the  transmitter  was 
the  same  as  used  in  the  homogeneous  case.  The 
results  for  this  case  are  presented  in  Figure  8  as 
propagation  loss  versus  range  at  a  receiver  height  of 
914  m.  The  overall  comparison  of  this  modeled  case 
and  the  measurements  is  quite  good  at  all  ranges. 
Similar  range-dependent  refractivity  structure, 
characterized  by  a  duct  that  rises  for  a  few  hundred 
kilometers  and  then  levels  out,  has  also  been  ob¬ 
served  in  the  San  Diego  area. 


Many  other  combinations  of  range-dependent 
ducting  could  be  modeled,  and  perhaps  a  fit  to  the 
measured  data  better  than  Figure  8  could  be 
achieved,  however,  since  definitive  meteorological 
measurements  wc'  not  made,  it  will  never  be 
possible  to  determine  with  certainty  what  propaga- 


RANGE  (km) 

Fig.  8.  Propagation  loss  versus  range  from  RK>  for  the 
elevated  duel  environment  of  May  28,  1974.  with  duct  rising  at  a 
slope  of  1:7,S0  to  2S0  km  and  no  change  beyond  250  km  {solid 
line)  compared  to  experimental  measuremcnis  (doited  line). 
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tion  mechanisms  were  responsible  for  the  relatively 
high  signal  levels  measured.  The  modeling  pre¬ 
sented  here  appears  to  confirm  the  conjecture  of 
Papperi  and  Good  hart  that  the  high  signals  resulted 
(rom  weakly  attenuated  modes  associated  with  the 
elevated  duct,  but  modified  by  mode-conversion,  or 
range -dependent,  effects.  These  results  also  con¬ 
firm  a  proposal  by  Wait  [1968]  that  the  usually 
neglected  whispering  gallery  modes  can  substan¬ 
tially  contribute  to  the  total  field  in  transhorizon 
propagation.  In  any  case,  the  modeling  presented 
here  seems  to  be  a  much  more  likely  explanation  of 
the  high  signals  than  the  evaporation  duct  model 
presented  by  Hitney  et  al.  [1978]. 

This  study  suggests  that  elevated  ducts  may  be 
more  important  than  often  considered  at  ranges  far 
beyond  the  horizon.  A  recent  investigation  by  Hit¬ 
ney  [1991]  used  RPO  in  a  parametric  study  of 
elevated  ducting  scenarios  to  infer  trapping  layer 
height  from  long-term  over-the-horizon  propagation 
measurements  at  547  MHz,  with  quite  good  results. 

CONCLUSIONS 

It  is  concluded  that  the  signal  levels  measured  on 
May  28,  1974,  at  the  greater  ranges  most  likely 
resulted  from  whispering  gallery  effects  in  an  ele¬ 
vated  duct  that  sloped  upward  with  range  away 
from  the  transmitter. 
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